Pathways to non-sequential multiple ionization in strong laser fields 
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The outgoing electrons in non-sequential multiple ionization in intense laser fields are strongly 
correlated. The correlations can be explained within a classical model for interacting electrons in 
the presence of the external field. Here we extend the previous analysis for two and three electrons 
to cases with up to eight electrons and identify the saddle configurations that guard the channels 
for non-sequential multiple ionization. For four and fewer electrons the electrons in the dominant 
configuration are equivalent, for six and more electrons this is no longer the case. The case of five 
electrons is marginal, with two almost degenerate transition configurations. The total number of 
configurations increases rapidly, from 2 configurations for three electrons up to 26 configurations for 
eight electrons. 
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I. INTRODUCTION 

Laser pulses with peak intensities of the order of 
10 14 W/cm 2 and wavelengths in the infrared can ion- 
ize several electrons simultaneously El 0, IdL Growing 
experimental evidence Jjj M , B p , 0, 0, El El El supports 
a rescattering process [La . fl3j as the most likely mech- 
anism for the rapid transfer of a sufficient amount of 
energy to the electrons. The interaction between the 
rescattered electron and the others takes place close to 
the nucleus, leading to the formation of a short lived, 
highly excited compound state with electrons near to the 
nucleus. This compound state will decay by emitting a 
single electron (the most likely decay mode) or several 
electrons. The highly excited electrons move so quickly 
that the changes in the field can be neglected: the elec- 
trons see an essentially static electric field. Experimen- 
tal studies of the final momentum distribution for double 
ionization show that the electron momenta parallel to the 
field are strongly correlated d, |tj El [H| . Building on 
the picture of an intermediate compound state we have 
proposed a classical origin for the electron correlations 
[UddllllLLi- If the energy of the compound state 
is above the zero field threshold energy no constraints 
on the outgoing dynamics are imposed. If the energy is 
below double ionization can occur only if the field is still 
present. Then the electrons will aim for the Stark saddle, 
but the electron repulsion will prevent them from cross- 
ing the saddle in the same place at the same time. Fur- 
thermore, if one electron is ahead it will push the others 
back, thus preventing multiple ionization. Therefore, the 
electrons have to escape in a configuration that balances 
the mutual repulsion. For two electrons this means that 
the electrons move side by side with respect to the field 
axis. Within this symmetric subspace it is easy to iden- 
tify a saddle which the electrons have to cross on their 
way to double ionization: it is a stationary configuration 
of two electrons side by side close to the Stark saddle for a 
single electron (assuming, as mentioned before, that the 



field is stationary). The number of saddle configurations 
increases with the number of electrons: there are two 
configurations for three electrons, with different energies 
and different threshold exponents [17|. 

We here extend our earlier study of such saddle con- 
figurations for two and three electrons to four or more. 
The aim is to identify the possible saddles, their critical 
energies and their threshold exponents. The number of 
configurations increases rapidly with the number of elec- 
trons. Beginning with five electrons also several nearby 
almost degenerate saddles appear. Most interestingly, for 
four and fewer electrons all electrons are equivalent in the 
dominant mode, but this is no longer true for six or more. 

We conclude this introduction with notation and the 
formulation of the Hamiltonian. Since the electrons es- 
cape from the highly excited compound state so quickly, 
we may assume that the field is constant during the es- 
cape. Let the electric field point in the z-direction and 
let the electrons be labeled i — 1, . . . , N with positions 
Yi = (xi,yi, Zi). Then the Hamiltonian is, in atomic 
units, 
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Simple scaling shows that the field strength F can be 
absorbed, so that without loss of generality we can take 

f = i EiGJ 

The aim now is to determine the stationary points of 
this Hamiltonian, where the derivatives with respect to 
momenta and positions vanish. For the most part this is 
not possible analytically. The numerical method of choice 
is the Newton- Raphson iteration jljj. The required cal- 
culation of the matrix of second derivatives can be turned 
to an advantage since it allows for an immediate deter- 
mination of the stability properties of configuration. The 
main disadvantage of the method, but also of all other 
numerical methods, is that they do not provide any mea- 
sure by which to judge whether all solutions have been 
found. On the practical side we can make sure that we 
have hit at least all the configurations with large domains 
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of attraction in the Newton-Raphson method by start- 
ing with a sufficiently large number of randomly selected 
initial conditions. We used up to 10 6 initial conditions. 
The fact that fairly similar but definitely different config- 
urations could be detected may be interpreted as strong 
evidence that all configurations have indeed been found. 

The stability analysis of the stationary point in the 
full phase space allows us to determine the behavior of 
the cross-section for non-sequential escape close to the 
threshold [Tsj . Among all Lyapunov exponents of the 
saddle there is one, X r , whose eigenvector components 
point in the same direction along the field axis: it corre- 
sponds to a simultaneous motion of all electrons in the 
same direction away from the saddle. Borrowing ter- 
minology from chemical reactions, we call this subspace 
the reaction coordinate. Because of the repulsion be- 
tween electrons all configurations have additional unsta- 
ble eigenspaces, which enter in the threshold exponent. 
If the initial energy of the system equals the saddle en- 
ergy only a trajectory living in the symmetric subspace 
can lead to a simultaneous escape of all electrons. This 
reduces the dimensionality of the problem and the cross- 
section vanishes. For higher energy some deviations from 
the symmetric motion are possible, giving a finite volume 
of initial conditions and a non-vanishing cross-section. 
The dependence of the cross-section on energy a(E) close 
to the saddle energy En can be obtained in the spirit of 
the Wannier analysis [2£j wl HE H3 , resulting in 



a(E) ~(E- E N y 



with an exponent 
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The Xi 's are the positive Lyapunov exponents of the sad- 
dle except for the reaction coordinate exponent A,.. The 
cross section is large if the exponent is small, i.e. if the 
saddle is crossed quickly (large A r ) or if the differences 
from the symmetric motion grow slowly (small Xi). These 
cross section exponents are an additional characteristic of 
the multiple ionization process. 

The remainder of the paper is organized as follows: 
in the next section we discuss a few obvious, and highly 
symmetric configurations. This is followed by our numer- 
ical results for up to eight electrons in section IIIII We 
conclude with a few remarks in section llVl 



motion to a dynamically allowed subspace in the high- 
dimensional TV-body phase space. That is, if in the full 
phase space of the TV-body problem initial conditions are 
prepared in this subspace they will never leave it. In the 
symmetry subspace we can use cylindrical coordinates for 
the electrons, Zi = z, pi = p and ipi = 2iri/N. Then, for 
zero total angular momentum along the field axis and in 
the scaled variables where F = 1, the Hamiltonian 
can be reduced to 



H(p P ,Pz,p,z) 
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The Hamiltonian I0J possesses a stationary point at p p 
0, p z = 0, p = pn and z = zn where 
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The first few saddles with electrons on a ring are shown 
in Fig.Hl 

For k <C TV the repulsion energy W behaves like 
W = (N 2 /tt) XX 1 A)- This suggests that W increases 
like TV 2 In TV, and indeed numerical evidence supports 



W w (0.37V 2 + 0.3iV-3.1)ln7V. 



(8) 



Since W increases faster than quadratically, the configu- 
rations exist for finite TV only. Specifically, we find that 
the energy of the configuration first decreases but then 
increases again, and that there is no configuration for 
N = 473 or more particles [35j. The properties of the 
configurations, including stability and cross section ex- 
ponents, are listed in table [I] 



One electron in the center 



II. SYMMETRIC CONFIGURATIONS 
A. All electrons on a ring 

We assume that all electrons are situated in a plane 
perpendicular to the field and that they obey a Cnv sym- 
metry. The reflection symmetry limits the momenta to 
be parallel to the symmetry planes and thus confines the 



The simplest extension from the symmetric ring is to 
place one electron on the field axis and the others on a 
ring as shown in Fig. Such configurations have three 
characteristic lengths, the radius of the ring pn-i, the 
position of the ring zat_i and the position of the electron 
on the axis z c . No analytical solutions have been found. 
The numerically determined characteristics are collected 
in table |n| 
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FIG. 1: Saddles with electrons on a ring. Shown are the cases 
with N from 3 up to 8, projected onto the x-y-plane perpen- 
dicular to the field. Note the increasing spacing between the 
electrons that indicates the growing influence of the electron 
repulsion. 



TABLE I: Saddles with electrons on a ring in a plane perpen- 
dicular to the field axis. The columns are the number N of 
electrons, the energy En at the saddle, the number n u of un- 
stable directions (excluding the one along which all electrons 
escape, i.e. the reaction coordinate), the Lyapunov exponent 
A r along the reaction coordinate, the critical exponent p near 
threshold, and the radius pN and the distance zn of the sad- 
dles from the nucleus. 



N 


En 


n u 


A r 




pN 


ZN 


2 


-4.5590 


1 


1.2139 


1.2918 


0.6580 


1.1398 


3 


-7.6673 


2 


1.1054 


2.6226 


0.9036 


1.2779 


4 


-11.1059 


3 


1.0340 


4.0971 


1.0994 


1.3882 


5 


-14.8004 


4 0.9817 


5.7342 


1.2677 


1.4800 


6 


-18.7044 


8 


0.9409 


8.6800 


1.4175 


1.5587 


7 


-22.7859 


10 


0.9077 


12.0849 


1.5538 


1.6276 


8 


-27.0208 


12 


0.8799 


15.9050 


1.6794 


1.6888 
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FIG. 2: Saddles with one electron on the field axis and the 
others on a ring. The electron in the center is shown shaded. 
Shown are cases N = 4, 5, 6, 7 and 8, projected onto the 
x-z-plane with the field axis along z (top row) and onto the 
x- y- plane perpendicular to the field (bottom row). 



TABLE II: Saddles with one electron on the field axis and 
the other N — 1 electrons on a ring. The columns are the 
number N of electrons, the energy En at the saddle, the 
number n u of unstable directions (excluding the one along 
the reaction coordinate), the Lyapunov exponent A r along the 
reaction coordinate, the critical exponent p near threshold, 
the position z c of the electron on the axis and the radius 
pjv-i and the distance zn-i from the nucleus for the ring of 
N — 1 electrons. 



N 


En 


Tin Aj- 




Zc 


pN-l 


ZN-l 


4 


-10.9398 


3 1.0271 


4.2423 


1.6543 


1.2996 


1.2718 


5 


-14.8001 


4 0.9750 


5.6633 


1.8199 


1.4324 


1.3950 


6 


-18.8975 


5 0.9345 


7.2035 


1.9690 


1.5583 


1.4960 


7 


-23.1867 


6 0.9015 


8.8984 


2.1054 


1.6775 


1.5813 


8 


-27.6372 


9 0.8739 


10.9698 


2.2317 


1.7905 


1.6553 



TABLE III: Saddles with all electrons in a plane through the 
field axis. The columns are the number N of electrons, the 
energy En at the saddle, the number n u of unstable direc- 
tions (excluding the one along the reaction coordinate), the 
Lyapunov exponent A r along the reaction coordinate, the crit- 
ical exponent p near threshold, and the minimal and maximal 
positions along the field axis, z m in and z max , respectively. 



N 



En n u 



3 -7.3902 3 1.0981 3.7040 1.1143 1.4666 

4 -10.3975 5 1.0177 6.8831 1.0295 1.5699 

5 -13.5302 7 0.9568 10.7953 0.9170 1.7755 

6 -16.7566 9 0.9082 15.4282 0.7898 1.8639 

7 -20.0549 11 0.8680 20.7764 0.6543 2.0199 

8 -23.4091 13 0.8339 26.8366 0.5141 2.0986 



C. All electrons on a line 

A third class of configurations has all electrons in a 
plane through the field axis. The electrons are placed 
like beads on a string. With the confinement of electrons 
to an almost linear arrangement the saddle energies are 
higher than in the previous cases and the cross section 
exponents are also much larger. The properties of the 
saddles are listed in table II I II and their positions are 
shown in Fig. [3] 



III. ALL CONFIGURATIONS 

A summary of the symmetric configurations discussed 
in the previous section is provided in Fig.^ which shows 
the saddle energies for different N. For the range of N 
included in the figure the saddle energies for all config- 
urations decrease, but as the example of the ring con- 
figurations shows, they can be expected to increase for 
sufficiently large N. There is a crossover in the config- 
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FIG. 3: Saddles with all electrons in a plane. Shown are 
cases N = 3, 4, 5, 6, 7 and 8, projected onto the :r-z-plane, 
with the field pointing along z axis. 
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FIG. 5: All saddles for four electrons, projected onto the x- 
z-plane with the field axis (top row) and onto the x-y-planc 
perpendicular to the field. Some electrons are shown shaded 
in order to help identify them in the two sets of frames. The 
states have Cav {y = 1), Cz v (v = 2), C v (y = 3) and Civ 
(y = 4) symmetries. 




FIG. 4: Energy of the saddles divided by N in the different 
symmetry classes for increasing number of electrons. The 
open circles are for electrons on a ring (Cnv symmetry), the 
crosses for configurations with one electron on the axis and the 
others on a ring (C(at-i)u symmetry) and full circles are for 
the configurations with all electrons on a line (C v symmetry). 



uration that gives the minimum, from the ones with all 
particles on a ring for N < 5 to the ones with a one 
electron in the center for N > 5. The case N = 5 is 
marginal (see below). In order to see what other config- 
urations one has to consider, we turn to a discussion of 
configurations found with the Newton-Raphson method 
for electron numbers between four and eight. The cases 
of two and three electrons have been analyzed previously 

0EEH0. 



TABLE IV: Saddles for four electrons, ordered by increas- 
ing saddle energy. The columns give the number v of the 
state, the energy E v at the saddle, the number n u of unstable 
directions (excluding the one along the reaction coordinate), 
the Lyapunov exponent A r along the reaction coordinate, the 
critical exponent [i near threshold, and some comment on the 
states. 



V 


E v 


n u 


A r /x comment 


1 


-11.1059 


3 


1.0340 4.0971 all on a ring 


2 


-10.9398 


3 


1.0271 4.2423 a ring plus center 


3 


-10.9398 


4 


1.0271 4.2708 - 


4 


-10.3975 


5 


1.0177 6.8831 all on a line 






A. 


Four electrons 



For four electrons there are 4 configurations, shown in 
Fig- El Their properties are listed in table IIVI The domi- 
nant non-sequential ionization takes place in the vicinity 
of the saddle corresponding to all electrons on a ring in 
a plane perpendicular to the field axis, as in the cases 
for two and three electrons. This saddle possesses both 
the lowest energy and the smallest critical exponent and 
thus the most favorable dependence of the cross-section 
on energy. 

The next two configurations are almost degenerate in 
energy and have similar critical exponents. They form 
a triangle with an additional electron placed in the cen- 
ter. The saddle v — 2 possesses a symmetry, but the 
state v = 3 has a C v symmetry only. The last v = 4 con- 
figuration, with all electrons in a plane, has significantly 
higher energy and larger critical exponent compared to 
the previous states. 
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FIG. 6: All saddles for five electrons, projected onto the x- 
z-plane with the field axis (top row) and onto the a;-?/-planc 
perpendicular to the field. Some electrons are shaded to help 
identification in the two sets of frames. The states possess 
Csv = 1), Cav = 2), C v (y = 3,4) and C 2v {v = 5) 
symmetries. 



TABLE V: Saddles for five electrons, ordered by increasing 
saddle energy. The columns give the number v of the state, 
the energy E v at the saddle, the number n u of unstable di- 
rections (excluding the one along the reaction coordinate), 
the Lyapunov exponent A r along the reaction coordinate, the 
critical exponent fi near threshold, and some comment on the 
states. 



V 


E v 




A* 


comment 


1 


-14.8004 


4 0.9817 


5.7342 


all on a ring 


2 


-14.8001 


4 0.9750 


5.6633 


ring plus center 


3 


-14.7763 


5 0.9766 


6.0823 




4 


-14.3922 


6 0.9684 


7.8093 




5 


-13.5302 


7 0.9568 


10.7953 


all on a line 



TABLE VI: Saddles for six electrons, ordered by increasing 
saddle energy. The columns give the number v of the state, 
the energy E v at the saddle, the number n u of unstable di- 
rections (excluding the one along the reaction coordinate), 
the Lyapunov exponent \ r along the reaction coordinate, the 
critical exponent [i near threshold, and some comment on the 
states. 



V 


Eu 




/" 


comment 


1 


-18.8975 


5 0.9345 


7.2035 


ring plus center 


2 


-18.7634 


6 0.9350 


7.7165 




3 


-18.7633 


7 0.9350 


7.7830 




4 


-18.7490 


6 0.9342 


7.9740 




5 


-18.7055 


7 0.9402 


8.5282 




6 


-18.7044 


8 0.9409 


8.6800 


all on a ring 


7 


-18.6476 


7 0.9293 


8.6653 




8 


-18.5511 


7 0.9297 


9.0837 




9 


-18.2526 


8 0.9241 


10.2688 




10 


-17.9228 


8 0.9205 


11.6055 




11 


-16.7566 


9 0.9082 


15.4282 all on a line 
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B. Five electrons 

For five electrons we find 5 configurations shown in 
Fig. H3 with properties given in table The five elec- 
tron problem constitutes actually a marginal case where 
the configuration of all electrons on a ring starts loos- 
ing its dominance. The energy of this configuration is 
only slightly lower than the energy of the state v = 2 
with one electron in the center and the others on a ring. 
The state v = 2 possesses, however, a smaller critical ex- 
ponent than the v = 1 configuration. Therefore, unless 
there is a wide disparity in prefactors in the cross section 
law or residual correlations in the compound state that 
prefer one state over the other, we expect that both will 
contribute significantly to the process of non-sequential 
five electron escape. 

The configurations v = 3 and v = 4 have the C v sym- 
metry and reveal higher energy and larger critical expo- 
nents than the previous ones. The last v = 5 configu- 
ration with all electrons on a line is the least significant 
because of its high energy and its high critical exponent. 



FIG. 7: The lowest lying saddles for six electrons, projected 
onto the x-z-plane with the field axis (top row) and onto the x- 
y-plane perpendicular to the field. Some electrons are shaded 
to help identification in the two sets of frames. The states 
have Csv (v = l), C v (v = 2), C 3v (y = 3) and C 2 v (v = 4, 5) 
symmetries. 



C. Six electrons 

For six electrons we find 11 saddles, with properties 
listed in table lVll The first five configurations are shown 
in Fig. The state with all particles on a ring is no longer 
dominant, it appears in the middle of table EU The 
saddle with an electron in the center and the others on 
a ring dominates the non-sequential ionization process. 
It has both the lowest energy and the smallest critical 
exponent. The states v = 2 and v — 3 look similar, 
but differ in their symmetries: state v = 2 has a C v 
symmetry, state v — 3 has a C^ v symmetry. 
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TABLE VII: Saddles for seven electrons, ordered by increas- 
ing saddle energy. The columns give the number v of the 
state, the energy E v at the saddle, the number n u of unstable 
directions (excluding the one along the reaction coordinate), 
the Lyapunov exponent X r along the reaction coordinate, the 
critical exponent fi near threshold, and some comment on the 
states. 



V 


E v 


n u 


A r 


/" 


comment 


1 


-23.1867 


6 


0.9015 


8.8984 


ring plus center 


2 


-23.0773 


7 


0.8980 


9.3618 


— 


3 


-22.9904 


8 


0.8991 


10.1259 




4 


-22.8409 


8 


0.8980 


10.7743 




5 


-22.8360 


8 


0.8967 


10.7509 




6 


-22.8126 


9 


0.9036 


11.3732 




7 


-22.7859 


10 


0.9077 


12.0849 


all on a ring 


8 


-22.7451 


9 


0.8932 


11.2920 




9 


-22.5550 


9 


0.8904 


12.3631 




10 


-22.4109 


9 


0.8905 


12.8292 




11 


-22.0806 


10 


0.8861 


13.8150 




12 


-22.0804 


11 


0.8861 


13.8850 




13 


-21.5215 


10 


0.8805 


16.1329 




14 


-20.0549 


11 


0.8680 


20.7764 all on a line 
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FIG. 8: The lowest lying saddles for seven electrons, pro- 
jected onto the 2-z-plane with the field axis (top row) and 
onto the z-y-plane perpendicular to the field. Some electrons 
are shaded to help identification in the two sets of frames. The 
states possess C§ v (y — 1) and C v (v = 2, 3, 4, 5) symmetries. 



D. Seven electrons 



For seven electrons we have 14 states, with the prop- 
erties listed in tabic IVlfl In Fig. Othe first five states are 
shown. The configuration with one electron in the cen- 
ter and the others on a ring dominates the simultaneous 
electron escape, just as for the six electron problem. 
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FIG. 9: The lowest lying saddles for eight electrons, pro- 
jected onto the :r-z-plane with the field axis (top row) and 
onto the a;-j/-plane perpendicular to the field. Some electrons 
are shaded to help identification in the two sets of frames. The 
states possess C% v {y = 1,2), C v (u = 3,4) and Crv (y = 5) 
symmetries. 

E. Eight electrons 

For eight electrons we have 26 configurations, listed in 
tablc lVlfll The dominant configuration has two electrons 
in the center and the other six around it, see Fig. EI The 
v = 2 saddle also has two electrons in the center, but 
different arrangements of the outer electrons. 

The state with one electron in the center and the others 
on a symmetric ring is the fifth configuration, the one 
with all electrons on a ring the 17th configuration. The 
saddles v = 3 and v = 4 are similar to the v — 5 one but 
they only have a C v symmetry, not the Ci v one as for 
state v — 5. All of these configurations possess similar 
characteristics, see table IVILTl 

The state with all particles on a line has the highest 
energy and the largest critical exponent. This is the case 
not only for the eight electron problem but for all N 
values analyzed here. 

IV. CONCLUDING REMARKS 

The analysis presented here confirms a trend antici- 
pated from our analysis of triple ionization: that with 
an increasing number of electrons more and more saddle 
configurations and thus channels for non-sequential mul- 
tiple ionization have to be considered. But the analysis 
also shows that, except for the somewhat marginal case 
of fivefold ionization, there is a well separated lowest ly- 
ing saddle that dominates non-sequential multiple ioniza- 
tion. Table iLXl collects the symmetries and the threshold 
exponents of this lowest lying state. 

A few experiments with multiple ionization of atoms 
have already been performed. However, momenta of the 
ionization products have been measured for double and 
triple ionizations only. Currently, measurements of mo- 
menta of more than two electrons are not available, not 
even of a single component, say the one along the field 
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TABLE VIII: Saddles for eight electrons, ordered by increas- 
ing saddle energy. The columns give an index v of the state, 
the energy E v at the saddle, the number n u of unstable di- 
rections (excluding the one along the reaction coordinate), 
the Lyapunov exponent A r along the reaction coordinate, the 
critical exponent fi near threshold, and some comment on the 
states. 



V 


E v 


n u 


Ay 




comment 


1 


-27.6592 


7 


0.8699 


10.6553 two in the center 


2 


-27.6523 


8 


0.8698 


10.9875 


_ 


3 


-27.6373 


7 


0.8737 


10.7629 





4 


-27.6373 


8 


0.8737 


10.7908 





5 


-27.6372 


9 


0.8739 


10.9698 


ring plus center 


6 


-27.6363 


8 


0.8724 


10.9400 




7 


-27.5495 


9 


0.8679 


11.5053 





8 


-27.5125 


8 


0.8683 


11.4266 





9 


-27.5122 


9 


0.8684 


11.5939 





10 


-27.4471 


9 


0.8667 


12.1438 





11 


-27.2884 


10 


0.8677 


12.9280 





12 


-27.2785 


10 


0.8692 


13.1026 





13 


-27.2537 


10 


0.8650 


13.4763 




14 


-27.2015 


10 


0.8665 


13.4391 




15 


-27.1801 


11 


0.8629 


13.4369 




16 


-27.0793 


11 


0.8741 


14.7803 




17 


-27.0208 


12 


0.8799 


15.9050 


all on a ring 


18 


-26.9869 


10 


0.8639 


14.3969 




19 


-26.9786 


10 


0.8652 


14.4260 




20 


-26.9774 


11 


0.8659 


14.7100 




21 


-26.9033 


11 


0.8608 


14.9062 




22 


-26.5081 


11 


0.8568 


16.7984 




23 


-26.3255 


11 


0.8567 


17.2700 




24 


-26.0390 


12 


0.8534 


18.2821 




25 


-25.1700 


12 


0.8463 


21.3654 




26 


-23.4091 


13 


0.8339 


26.8366 


all on a line 



preference for similar momenta. For N — 5 and more 
escaping electrons additional structures in the momen- 
tum distribution should appear, since the lowest lying 
saddle has non-equivalent electrons, see Figs.|51to[!2] The 
momentum distributions for the two N — 3 saddles are 
discussed in |T^ . and similar modifications in shape can 
be expected for more electrons. 

The analysis presented here assumes simultaneous es- 
cape of all electrons in a multiple ionization event. It is 
of course also possible that first one or a few electrons 

TABLE IX: Saddles with lowest energy for different number 
of electrons. The columns give the number TV of electrons, 
the energy En at the saddle, the critical exponent n near 
threshold, the total number of states vn and some comment 
on the states. 



N 


En 




vn comment 


2 


-4.5590 


1.2918 


1 all on a ring 


3 


-7.6673 


2.6226 


2 all on a ring 


4 


-11.1059 


4.0971 


4 all on a ring 


5 


-14.8004 


5.7342 


5 all on a ring 


6 


-18.8975 


7.2035 


11 ring plus center 


7 


-23.1867 


8.8984 


14 ring plus center 


8 


-27.6592 


10.6553 


26 two in the center 



escape, then the next batch of electrons leaves and so 
on, before one arrives at the final multiply charge ion. 
Sequential ionization becomes significant when the field 
intensity is high enough to ionize electrons from a highly 
excited state of the ion left behind when the first elec- 
trons are gone. Nevertheless, close to the threshold for 
non-sequential multiple ionization the characteristics of 
the non-sequential process as described here should be 
identifiable. 



axis. But already the measurement of parallel momenta 
for just two electrons would suffice for a preliminary test 
of the ionization pathways considered here. In particular, 
the breaking of symmetry between the outgoing electrons 
should be accessible. As long as the electrons are inter- 
changeable, the distribution of parallel electron momenta 
should be like that for double ionization, with a clear 
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